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ABSTRACT 

We present rest-frame 15 and 24 /xm luminosity functions and the correspond- 
ing star-forming luminosity functions at z<0.3 derived from the 5MUSES sample. 
Spectroscopic redshifts have been obtained for ~98% of the objects and the me- 
dian redshift is ~0.12. The 5-35 /xm IRS spectra allow us to estimate accurately 
the luminosities and build the luminosity functions. Using a combination of 
starburst and quasar templates, we quantify the star-formation and AGN con- 
tributions in the mid-IR SED. We then compute the star-formation luminosity 
functions at 15 fim and 24 /im, and compare with the total 15 /xm and 24 /xm lu- 
minosity functions. When we remove the contribution of AGN, the bright end of 
the luminosity function exhibits a strong decline, consistent with the exponential 
cutoff of a Schechter function. Integrating the differential luminosity function, 
we find that the fractional contribution by star formation to the energy density is 
58% at 15 /im and 78% at 24 /im, while it goes up to ~86% when we extrapolate 
our mid-IR results to the total IR luminosity density. We confirm that the ac- 
tive galactic nuclei play more important roles energetically at high luminosities. 
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Finally, we compare our results with work at z~0.7 and confirm that evolution 
on both luminosity and density is required to explain the difference in the LFs 
at different redshifts. 

Subject headings: galaxies: active - galaxies: starburst - galaxies: luminosity 
function - infrared: galaxies 



1. Introduction 



The unprecedented sensitivity of the Spitzer Space Telescope (IWerner et al.l 120041 ) has 
opened a new window to e xplore the infrared (IR) unive rse. IR lurninous galaxies, dis- 
covered from ground-based (IRieke fc Lowlll972l ) and space (ISoifer et al.lll987l ) observations, 
constitute an important population. While UltraLuminous InfraRed Galaxies (ULIRGs, 
LiR > IO^^Lq) only account for ~5% of the IR luminosity in the local universe, their contri- 



bution becomes increasingly important at higher redshift, e.g. Luminous InfraRed Galaxies 
(LIRGs, lO^^Lp^ < Ltr < IO^^Lq) are responsible for 70%±15% of the energy density at z~l 



( iLe Floc'h et al.ll2005l ) and ULIRGs become more important than LIRGs at z~2. 



Recent surveys with the Multiband Imaging Photometer for Spitzer (MIPS) ( iRieke et al. 



20041 ). as well as earlier observations with the Infrared Space Observatory (ISO), probe the 
dust emission at much fainter levels than that has been reached by the Infrared Astro- 
nomical Satellite (IRAS). This has vastly improved our understanding of galaxy evolution. 
Surveys with I RAS first established the local benchmark for mid- and far- IR luminosity 



functions (LPs) (jSoifer et al.lll987l : ISaunders et al.lll990l : iRush et al.lll993l : IShupe et al.lll998l . 

Sanders et al. joosh. In the nineties. 15/im ISOCAM observations were used to derive the 



15 /im LP (IXul I2OOOI) at low redshift, as well as to s tudy the evolution e ffects from the 



numb er counts (lElbaz et al.lll999l : IChary fc Elbaal200ll ). The recent work of iBothwell et al. 



(120111 ) has constrained the slope of the IR and UV LPs at the extreme faint end for the first 
time using large datasets of local galaxies, and has derived the distribution function of star 
formation rate in the local universe. Deep MIPS surveys c arried out in the past fe w years fer- 
tilized the ground for understanding the evolution of LP. iLe Ploc'h et al.l (120051 ) illustrated 
the variation of 15 yum LP in the range of 0.3< z <1.2 and suggested that the comoving IR 
energy density evolves dramatically, in c reasing with look-ba ck time as (l+z)^'^^"'^ up to z~l 
JCaputi et al.lboO?! : iReddv et allboosl : iMaenelh et aDl2009h . 



IR bright galaxies emit the bulk of their energy as dust-reprocessed light generated by 
dusty star formation (SP) or accretion onto the supermassive blackholes, referred to here- 
after as active galactic nuclei (AGN). Obtaining information on the relative contribution of 
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SF/AGN is critical for understanding a galaxy's integrated emission. iLe Floc'h et al.l (|2005[ ) 
explored the star-formation history at 0.3<z<1.2. However, the MIPS 24 /xm photometry 
they used for that study did not allow them to account for the AGN contribution, or con- 



(2010 


); 


Fu et al. 


(2010): 


Goto et al. 


(2011) 



have identified AGN dominated sources and excluded these objects from their samples to 
derive the SF LFs. The shape of LF depends on th e rest-frame wavelength. UV/optical LF 



usually follows the profile of a Schechter function (lArnouts et al.l l2005l : lllbert et al.ll2005l ). 



while at mid/far-IR wavelengths, the brig ht end slope is observed to be shal l ower than the 



expo nential cut-off of a Schechter function (jSoifer et al.lll987l : iRush et al.lll993l : ISanders et al. 



20031 ) . This is rather intriguing because both UV and far-IR emission traces the same star- 
formation, and the differe nt shapes for the LFs are suggested to be a result of dust extinction 
effect. The recent work of iFu et al.l (120101 ) claims that the shallower slope of the IR LF could 
be due to the contribution of AGN at the high luminosity end and when the AGN component 
is removed, the SF LF can be fit with a Schechter function. This further motivates our work 
of separating the SF and AGN emission in our objects to understand their contribution to 
the LF. 

To quantify the relative contribution of star-formation and AGN to the infrared lumi- 



nosities, spectral decompositions have been performed by several group s ( ISajina et al.ll2008 



Pope et al.ll2008l : Murphy et al.l 120091: iMenendez-Delmestre et al.ll2009l ). In high mid-IR lu- 
minosity systems at z~2, ISajina et al.l ( 120081 ) found an average AGN fraction of ~ 20% — 30% 
of total IR luminosity for st rong PAH sources while this number goes up to >70% for weak- 
PAH sources. The study of iPope et al.l (120081 ) of a sample of submm galaxies at similar red- 
shift has revealed at most a 30% contribution from AGN at mid-IR (5-11.5 /im rest-frame). 



In the local universe, the contribution of nuclear activity to the bolometric 



umin osity of 
(l2009l) . reach- 
J20081 ) suggest 



ULIRGs has been quantified with six independent methods b v lVeilleux et a 
ing an average AGN contribution of ~ 35% — 40%, whereas iNardini et al. 
that intense star formation accounts for 85% of the IR emission in local ULIRGs, with AGN 
contributing 15%. It i s clear that the relative contribution o f SF/AGN varies in galaxies of 
different luminosities (lYuan et al.ll2010l : iHopkins et al.l 120101 ) . with the AGN playing more 
critical roles in more luminous systems. It must also be a function of wavelength, since the 
spectral energy distribution (SED) of the two components are quite distinct. The different 
selection criteria and decomposition techniques among different samples and authors also 
add to the differences derived in the relative contribution of SF and AGN. 

Current studies on IR LFs are mostly based on MIPS 24 /im observations, which implies 
a heavy reliance on the SED library used to make k-c orrections and d erive the monochromatic 
continuum or total infrared luminosities. Recently, iFu et al.l (120101 ) have derived 8 fim and 
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15 yum LFs, as well as the corresponding star-formation LF, for a sample of z~0.7 objects, 
taking advantage of their IRS spectra. The 5 Millijansky Unbiased Extragalactic Survey 
(5MUSES), an infrared spectroscopic survey of 330 objects selected by their MIPS 24 /xm 
flux densities, provides an important sample for understanding the infrared galaxy population 
(Helou et al. 2011, in preparation). Although IR luminosity functions have been extensively 
studied at low redshift, we find 5MUSES to be a unique sample for deriving LF for the 
following reasons: 1. Being a mid-IR flux-limited sample, we reach a wide range of infrared 
luminosity (^lO^L^ to 2. We reach to higher redshifts than previous samples, e.g. 



the RBGS (ISanders et al.ll2003l ). 3. The IRS spectra of our sample and far-IR measurements 
for most of them allow us to minimize the uncertainties on k-correction and luminosity 
estimation associated with adopting specific SEDs to be applied to all sources. 4. Last 
but not least, the IRS spectra allow for a careful decomposition of every single source into 
a SF and AGN component and estimation of their contribution to the luminosity density. 
This will facilitate further studies on how the SF/AGN fraction in the integrated luminosity 
density evolves in a cosmological context. 

The paper is organized as follows. In Section 2, we briefly describe the sample selection, 
and then the infrared and optical data used in this study. We introduce our methodology of 
using the 1/Vmax method to derive luminosity function in section 3, where we also demon- 
strate how we correct for the incompleteness of the 5MUSES sample. In the same section, 
we discuss in detail our methods of the SF/AGN decomposition of the IRS spectra and how 
we estimate the SF contribution in a statistical sense. The 15 and 24 fim LFs, as well as the 
corresponding SF LFs are presented in Section 4, together with a discussion of how the star- 
formation fraction varies with luminosity and wavelength. We summarize our conclusions 
in Section 5. Throughout this paper, we assume a ACDM cosmology with Hq =70kms~^, 
Q^=0.27 and Qx=0.73. 



2. Observation and Data Reduction 



2.1. The Sample 



5MUSES is a mid-infrared spectroscopic survey of 330 galaxies with 24 fim flux densities 
5mJy<fy(24/im)<100 mJy. The sources are selected from the SWIRE (Elais-Nl, Elais-N2, 
Lockman Hole and XMM) and the Extragalactic First Look Survey (XFLS) fields, covering 
a total area of 40.6 square degrees on the sky. It provides a representative sample at interme- 
diate redshift {{z) = 0.144), previously unexplored by Svitzer since most of the spectroscopic 
work was foc used on nearby spiral galaxies (SINGS) (IKennicutt et al.ll2003l ). local LIRGs 



and ULIRGs flArmus et al.l 120071 . 120091 : 1 Veilleux et al.ll2009[ ). and the much fainter and more 
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distant (z~2) galaxies ( iHouck et al.l 120051 : lYan et al.l 120071 ). A total of 1111 objects have 
fi,(24/im) between 5 and 100 mJy from the five survey fields of 5MUSES excluding stars. In 
order to efficiently observe the objects using the staring mode of IRS and include the largest 
fraction of a galaxy's integrated light, only objects unresolved within an aperture of d=10.5" 
(corresponding approximately to the slit width of the Long-Low module of IRS) are included 
in the final pool and this results in a total of 800 sources. Then 330 objects are randomly 
selected from the 800 final candidates. The details of the sample can be found in Helou et 
al. (2011, in preparation). 



2.2. Data Reduction 

We have obtained low-resolution mid-IR spectra for all 330 objects in 5MUSES using 
the Infrared Spectrograph on board the Spitzer Space Telescope. Both Short-Low (SL: 
5.2-14.5 /xm) and Long-Low (LL: 14-37/im) modules are used, with spectral resolution of 
64-128. The integratio n time on each object ranges from 300 to 960 seconds to achieve 



roughly the same SNR ( iWu et al.ll2010f ). The IRS data are processed by the Spitzer Science 
Center data reduction pipeline version S17 and our data reduction starts from the pipeline 
products designated as "basic calibrated data (bed)". The two-dimensional spectrograms are 
median combined and then the off-source sky regions are subtracted. After removing the sky 
background, the spectrograms are cleaned with the IRSCLEAN package to remove bad pixels 
and apply rogue pixel correction. Then the background-subtracted cleaned spectrograms are 
reduced with the optimal extraction method of the Spitzer IRS Custom Extractor (SPICE) 
software to extract the 1 -D spectra. The details on the reduction of the IRS data can be 



found in IWu et al.l fcOlOf ) . 



The IRS spectra allow us to derive redshifts for sources with mid-IR emission and/or 
absorption features. We have also searched for optical spectra of our sample in the literature. 
A total of 50 5MUSES objects either do not have redshift information from the literature, 
or their IRS spectra are characterized by featureless power law continua, which cannot yield 
redshifts. We obtained redshifts for some of these sources by using the double spectrograph 
instrument on the Palomar 200 inch Telescope in 2009B and 2010A terms. The optical 
data are reduced with the IRAF software using the standard routines for bias subtraction, 
flat-fielding, sky background removal and wavelength calibration. Finally, for sources with 
multiple emission line features, redshifts are derived using all available emission lines. 
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2.3. The Redshift Completeness 

Combining the IRS, Palomar spectra and literature work, we found secure redshift mea- 
surements for 309 objects out of the 330 sources in 5MUSES (94%). The redshift distribution 
of 5MUSES is shown in Figured] We also show the redshift distribution for starburst, com- 
posite and AGN-dominated sources respectively, which have been classified based on their 
apparent 6.2 /im PAH EWsQ (Wu et al. 2010). Since we are only interested in the red- 
shift range of z<0.3 for this study, the relevant redshift completeness is close to 1 for two 
reasons: (1) The sources for which we are not able to find redshifts from 80 minutes of 
integration time on the Palomar 200 inch telescope have very low r band to 24 /im band 
flux ratios fi,(r)/fi,(24/im). As can be seen from Figure |21 objects with low fj,(r)/fj,(24/im) 
ratios (ie. log[f,y (r)/f^(24jUm)]<-2.6) are more likely to have high redshifts (z>0.3) (see also 



Dey et al.ll2008l ): (2) The IRS spectra of the sources without redshifts are characterized by 
featureless power-law continuum in the mid-IR. This indicates that they are most likely 
AGN-dominated. It can be seen from Figured] that the median redshift for AGN-dominated 
sources is much higher than SB or composite sources. The median redshift for sources 
with 6.2 /im PAH EW> 0.2 /im is 0.13, while it is 0.40 for sources with 6.2 /im PAH EW < 
0.2 /im. Thus these power law sources are much more likely to lie at the high end of the 
redshift distribution. 

For this study, we focus on objects with z<0.3, which includes 226 objects. Among the 
21 (330-309) sources without redshift, only 4 do not have very low f^{T)/{^{24fim) ratios, 
e.g. their log[f,y(r)/f,^(24/im)]>-2.6, and might be located at z<0.3. This indicates that the 
redshift completeness for our sample at z<0.3 is >~98%. 



3. Methodology 

3.1. The Incompleteness Correction 

The targets for 5MUSES are randomly selected based on fi,(24/im)>5mJy after exclud- 
ing the resolved objects. Understanding the selection function for 5MUSES is crucial for 
building the luminosity function. Because we exclude extended objects, it is likely that we 
have excluded more nearby objects than those at higher redshifts. Thus when we derive the 
number density, instead of applying a uniform correction factor of ~3.4 (1111/330), we need 
to investigate the selection effect in individual redshift bins before we build our luminosity 



^Sources with 6.2 /im PAH EW>0.5/im are SB-dominated; sources with 0.2 /im<PAH EW<0.5/im are 
composite; and sources with PAH EW<0.2 fim are AGN-dominated. 
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function. 



Redshift information is not a vailable for all the sour ces in the parent sample of 5MUSES, 
thus we use the redshift catalog of lPapovich et al.l (120061 ) for the XFLS field to unde rstand the 
redshift distribution when a flux limit of 5 mJy is imposed. iPapovich et al.l (120061 ) observed 
the XFLS field using the Hectospec instrument on MMT for 5 positions covering a 1 degree 
diameter field of view individually. Then they combined redshifts from Hectospec with 
redshifts from SDSS, and reached a completeness of ~90% at fi,(24/xm)>l mJy in the 3.3 
square degrees of their survey field. These authors also provided the completeness factors at 
different flux limits, which we use to derive the final number counts at fj,(24/im)>5mJy in 
the XFLS field. Then we divide the number of objects in different redshift bins by the total 
number of objects in the XFLS field, and derive the fractional contribution of number counts 
in this field at fj,(24/im)>5mJy. Using this as a reference, we predict the number of objects 
in the corresponding redshift bins for the 5MUSES sample. Then we divide the predicted 
number counts by the number of objects we have observed and obtain the correction factor 
in each redshift bin. Finally, we fit a second-order polynomial to the data and this gives us 
the completeness correction factor Ui{z), which is then used to correct for the incompleteness 
at different redshifts. 



3.2. The 1/Vmax Method 

We limit our study of the mid-IR luminosity function to z<0.3 because the rest-frame 
24 /xm bancj^ moves outside the IRS spectrum beyond z=0.3. In addition to that, our rela- 
tively bright flux limit of f,y(24/im)>5mJy results in a fast decrease of the number of objects 
as redshift increases, which would yield results that have low statistical significance at high 
redshift. 



In this study, we use the 1/Vmax method (ISchmidtlll968l ) to derive the luminosity func- 



tion, which does not require any assumption on the shape of the LF. The l/Vmax method 
counts galaxies within a volume. Vmax is calculated individually for each source in our 
sample as the maximal volume within which that galaxy is detectable in this survey. The 
availability of 5-35 /xm IRS spectrum allows us to accurately make k-corrections based on the 
observed SED shape for individual galaxy. We first derive the k-correction for each object, 
and then move the galaxy to the redshift where its 24 fim observed flux reaches the limit of 
this sample, 5 mJy. The maximum comoving volume is calculated as: 



^Here we refer to the MIPS 24 /im filter, instead of the monochromatic 24 /im continuum. 
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"Zhigh 

Vi.max = / -rdz (1) 



OW 



dz 



where [ziow,Zhigh] is the redshift range of interest. For our study, we ignore sources with 
z<0.02, while 

Zhigh is the lower of the two: a) the maximum redshift considered in this study, 
0.3, b) the maximum redshift a source can be detected at the limiting observed 24 //m flux 
of 5mJy. The uniformity of the distribution of galaxies is tested by checking the V/Vmax 
values and we find (l^/Knax)=0.54 for the sources used in this study. 



The luminosity function is then derived by using the following formula ( 1Schmidtlll968l ): 

Where fl is the total survey area of 5MUSES sample (40.6 square degrees), Vi^max is the 
comoving volume over which the ith galaxy could be observed, AlogL is the size of the 
luminosity bin and Ui is the completeness correction factor for the ith galaxy. Ui is a function 
of redshift and was calculated in Section 3.1. We divide the sources into seven luminosity bins, 
and calculate the value of (j) in each bin. The uncertainties include both the Poisson noise 
statistics on the number of sources used in the measurement, and the uncertainty associated 
with the completeness correction factor cuj. As can be seen from Figure [3l the uncertainty of 
Ui is rather large, mainly due to the small number of objects in each redshift bin in XFLS at 
fiy(24/im)>5 mJy, so we assign an uncertainty of 40% to Ui, which is the average uncertainty 
for the data points we use to calculate the correction factor. Because our k-corrections are 
made directly from the source SED, we have almost negligible uncertainties associated with 
the conversion from the observed flux to the rest-frame luminosities. For the luminosity 
functions presented in this paper, we do not include noise from the cosmic variance since we 
sample several widely distributed directions. 



3.3. Spectral decomposition in the mid-IR 

The mid-IR SED of star forming galaxies and AGN show distinctly different spectral 
features. As a result, the availability of 5-35 fim IRS spectra for the 5MUSES sample al- 
lows us to disentangle the star-formation and AGN contribution to the galaxy luminosity. 
Star-forming galaxies often display broad emission features, which are generally attributed 
to the emission from PAHs. AGN, on the other hand, are usually characterized by feature- 
less power-law continuum (except for a few high-ionization fine-structure lines) and their 
mid-to-far IR continuum slopes are normally flatter than star-forming galaxies. A combina- 
tion of one starburst template and a power-law continuum with a free spectral index have 
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often been used to dec ompose galaxy sp e ctra into star-formation and AGN components for 



high-redshift galaxie s (ISajina et al.ll2008l : iPope et al.ll2008l : iMenendez-Delmestre et al.ll2009 



Murphy et al.l l2009f ). The high S/N IRS data of the 5MUSES sample allow us to take 
into account detailed mid-IR features and decompose the galaxy spectra much more accu- 
rately. Among mid-IR spectral features (PAH strength, continuum slopes, silicate strength, 
fine-structure lines, etc.), PAH strength or the IR continuum slopes are arguably the most 
commonly used indicators for star-formation, so we select our templates mainly based o n 



these two parameters. We start from our empirical SED template library (jWu et al.ll2010l ). 
which has included star-forming galaxies, ULIRGs and PG/2MASS quasars, and select 15 
star-forming galaxy (6.2 fim PAH EW>0.5 fim) templates and 10 quasar templates, cover- 
ing as large a range of PAH strength and slope variation as possible. Then we perform a 
least-x^ fit for each combination of a star-formation template and a quasar template to find 
the most likely coefficients that would describe the observed 5MUSES spectrum as a linear 
combination of the two templates. In the upper panels of Figure HJ we show examples of de- 
composition of IRS spectra of typical starburst dominated, composite and AGN dominated 
sources. 

Statistical constraints on the star-formation contribution: Rather than directly 
adopting the star- formation fraction at 15 fim from the best-fit, we build the probability 
density function {pdf) of the star formation fraction by weighting the values of star formation 
fraction for each trial fit by exp(-x^/2). Then the star-formation fraction is taken to be the 
median of the resulting probability density function and the la uncertainty is taken to be 
the 16th-84th percentile range. On the lower panels of Figure IH we show the corresponding 
pdf for each source and the SF fraction and its associated uncertainties are also indicated 
on the plot. As discussed earlier, the 6.2 fim PAH EW and the continuum slope have often 
been used as indicators of star-formation activities. We compare the 6.2 /xm PAH EW and 
continuum slope of fi,(24/im)/f,^(15/im) versus SF fraction estimated from the probability 
distribution on the left and middle panels of Figure O Our median likelihood SF fraction 
clearly correlates with both parameters while the scatter is quite significant for intermediate 
values of f^(24yum)/f^(15yum) ratios. Note however that for f,^(24/xm)/fi,(15/im)<0.3, all but 
2 sources have very low star-formation fraction, suggesting that very fiat mid-IR slopes are a 
strong discriminator for AGN. We also note that some sources, even though they have very 
large PAH EWs (>0.5 /um), have star-formation fractions only >~0.5. This could be due to 
the fact that we are looking at the star-formation fraction at 15 fim in this study, while the 
contribution to the total IR luminosity is very likely to be dominated by SF. On the other 
hand, we also need to point out that the SF fraction we derive has at least ~20% uncertainty. 
Finally, we plot the SF fraction versus the IR luminosity for each source on the right panel 
of Figure [5] and we do not observe any correlation. This indicates luminosity itself provides 
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limited information on the energy source of a galaxy. Clearly, there is no single parameter 
that could be used to determine the SF fraction accurately, thus a combination of a few is 
indeed needed in order to constrain the relative contributions of SF and AGN. 



Results 



4.1. The 24 /im and 15 /xm Luminosity Functions at z<0.3 



From the 5-35 /im IRS spectrum of the 5MUSES sample, we can directly estimate the 
rest-frame MlPS-equivalent 24 /xm luminosities. Using the 1/Vmax method, we have derived 
the 24 /im luminosity function for sources with z<0.3. Our 24 /im LF is shown in Figure [6] 
and the correspondin g data points are rep orted in Table [H We adopt a double power-law 
exponential function (ISaunders et al.lll990l ) to fit our LF: 



dN{L) 



dVdlogiQ{L) 



,L_ 



exp{-^logUl + {^)]} 



(3) 



The dashed line is the fit to our data by running mpfit.pro. The uncertainties on the 
fitting parameters are derived with 1000 Monte Carlo realizations and they are reported 
in Table [21 The dotted line denotes the luminosity calculated at the median redshift of 
our sample corresponding to the 24 /im flux limit. Note that 5mJy is the flux limit on the 
observed 24 /im, while the flux limit we use here is taken to be the rest frame 24 /im flux 
density corresponding to the galaxy with fnhs(24 Ai ni)=5 mJy and ma x imurn k-corrected. On 
Figure we have also included t he 24 /im LFs by iRuiopakarn et al.l (120101 ) for comparison. 
Our LF is in good agreement with Rujopakarn et al. (j201ol ) at the low end, while at the high 
end, our data points are located between the 0.05<z<0.2 and the 0.2<z<0.35 LFs they have 
derived. This indicates that in the redshift range we derive our LF, evolution is probably 
already at work, as will be addressed in more detail later on. 

We then compute the rest-frame 15 /im monochromaticEl luminosity function for our 
e. LF at this wavelength has been extensively studied in t he local universe (|Xu 



samp 
2000l: 



Pozzi et al.ll2004J : iMatute et al.ll2006f ) as well as at high redshift flLe Floc'h et al.ll2005l : 



Magnelli et al.ll2009l ) to explore the evolution of galaxy populations. We repeat the steps 
used for calculating the 24 /im LF and derive the 15 fim LF for the 5MUSES sample. In the 
upper panel of Figure [TJ we plot our 15 /im LF with solid circles and the corresponding data 
points are reported in Table [H The dotted line denotes the luminosity calculated at the 



■^Thc monochromatic luminosity is calculated within a width of 1 /im. 
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median redshift of this sample corresponding to the 15 fim flux limit. As expected, the com- 
pleteness level at 15 fim is much lower than at 24 /im, since this is a 24 fim selected sample 
at z<0.3. The dashed line is the fit to our LF with double power-law exponential function. 
During our fit, we fixed the faint end slope to q;=1.2, which has been well determined from 
similar studies in the local universe. The results of the fitting parametersare reported in 
Table [2l For comparison, we also overplot the local 15 /zm LF from IXul (120001 ) derived using 
ISO CAM observations. Simil ar to what we have observed when we compare our 24 /xm LF 
with iRujopakarn et al.l (120101 ). the 5MUSES LF is in good agreement with Xu's LF at the 
low end, while our LF is slightly higher at the high luminosity end. This can be explained by 
the differential evol ution effect. The density or luminosity evolution with redshift have bee n 
extensively studied (ILe Floc'h et al.ll2005l : IPerez- Gonzalez et al.ll2005l : iMagnelli et al.l 120091 ). 
Evolution affects different luminosity bins by different amounts, because they are populated 
from different redshift ranges. In the lower panel of Figure [TJ we plot the median redshift in 
each luminosity bin as a solid circle and the 16th and 84th percentile as the error bar. The 
median redshift (z=0.09) for Xu's sample has also been plotted (solid line) together with its 
Icr population dispersion. Although the median redshift for the sources used in this study 
is 0.12, it is clear that the high luminosity end is dominated by sources at higher redshift, 
thus our slightly higher LF at the bright end is almost certainly a result of that difference. 

In addition to using the double power-law to fit the LF, we have also attempted using 
a single Schechter function. In Figure U\ we overplot the fit with Schechter function as the 
blue dot-dashed line. Although Schechter function can fit our data reasonably well, the fit 
is rath er poor in both the faint and bright end when the data points from the ISO 15 fim 
LF byH J2000h are included. 



4.2. The Star-formation Luminosity Function 



In the previous subsection, we have shown that the 15 fim LF need to be fit with a double 
power-law exponential function because the bright-end slope of the LF is clearly flatter than 
the Schechter Function. This is in contrast to UV LFs, which display much steeper slopes 
at the high luminosity end. As both the UV and IR luminosities trace active star formation, 
with the IR being the reprocessed portion of the UV, one would expect similar behaviors 
by the LFs of UV and IR emission. Using Spitzer IRAC observations of the Bootes field, 
Huang et al.l (120071 ) have shown that the 8 fim luminosity function of a sample of star-forming 
galax ies does ind e ed fol low the shape of a Schechter function. More recently, using AKARI 
data, iGoto et al.l (|201ll ) have shown that after removing t he optical l y iden tified AGN, their 
IR LF becomes much steeper. Using Spitzer IRS spectra, IFu et al.l (|2010[ ) have studied the 
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15 yum LF at z=0.7 and constrained the slope of the SF LF at the high luminosity end using a 
Schechter function. So the flatter LFs we are deriving probably reflect AGN contributions to 
the IR, and we will verify this hypothesis in what follows. In Section 3.3, we have shown our 
method of decomposing the star- formation/ AGN contribution from the mid-IR spectra of 
5MUSES sources and estimated the star-formation fraction by taking the median likelihood 
of the probability density function. Now we derive the star-formation luminosity functions 
at 15 /im and 24 //m for 5MUSES at z<0.3. 

Local 15 /im SF LF: Using the median estimate of the star-formation fraction for each 
object, we derive the 15 /xm star-formation luminosity by multiplying the 15 /xm monochro- 
matic continuum luminosity by the star-formation fraction at the corresponding wavelength 
obtained from spectral decomposition. The resulting individual SF luminosities are then 
used to build the SF LF. On the left panel of Figure [HI we show the 15 /im SF LF (red 
squares) at z<0.3. The black dotted line is a fit to the LF adopting the Schechter function. 



dN{L) 



dVdlogio{L) 



(4) 



We overplot the local 15 fim LF for normal spiral and starburst galaxies by lPozzi et al. 



(120041 ) as the blue crosses. These authors analyzed data from the ELAIS southern fields and 
excluded AGN in their study. Their SFLF appear to be in good agreement with our results. 
We also overplot the total 15 /im LF as the black filled circles on the left panel of Figure 
[HI and the dashed line is a fit to the total 15 fim LF with double power-law exponential 
function. The two LFs, the total and the star-formation LFs at 15 /zm, differ significantly at 
Lis^jm > IO^^Lq. This is presumably due to the AGN contribution at high lum inosities. The 



AGN L Fs display disti n ctly d ifferent shapes from SF LF as has been shown by lMatute et al. 



(I2OO6I ): [Hopkins et al.l (120071 ). etc. When the luminosity increases, the AGN contribution 
also increases progressively. This AGN component reveals its presence in the total 15 /im LF 
by requiring a different slope in the fit at high luminosity. In Figure [HI we have overplotted 
the 15 /im AGN LF as green diamonds and fit the data with a double power-law exponential 
function. Clearly, the AGN LF presents a much shallower slope at high luminosities. On 
the other hand, the 15 /im SF LF drops quickly at the bright end and it could be fit well 
with a Schechter function. This has alre a dy be en seen in the local universe for the 8 /im LF 
of star formin g galaxies bv [Huang et al.l (|2007[ ). as well as the 15 /im AGN-corrected LF at 
0.6<z<0.8 bv lFu et all tom . 



Local 24 /im SF LF: We repeat the same steps used to derive the 15 /im SF LF at 
24 /xm and show the 24 /im SF LF on the right panel of Figure [HI Again, we observe the 
departure of the tot al and SF LFs at 2 4 /im a t the bright end and the best-fit parameters are 
reported in Table [2l [Rujopakarn et al.[ (120 10[ ) have identified AGN from optical spectroscopy 
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and derived 24 /xm SF LF in several redshift bins. Their comparison of the 24 fim total LF 
and the 24 fim SF LF shows a very similar trend as what we have observed from our sample. 
We also overplot their 24 yum SF LFs at 0.05<z<0.2 and 0.2<z<0.35 as the blue diamonds 
and green crosses on Figure [Si 



4.3. Discussion 

As the main luminous phenomena in the universe, star formation and AGN activities 
have been extensively studied across all wavelengths. Despite the amount of effort to quan- 
tify the SF/AGN contribution and explore its evolution with luminosity /redshift, no fair 
comparison can be made unless truly equivalent samples are studied. Our 5MUSES sample, 
after correcting for selection effects, essentially defines a relatively bright IR selected unbi- 
ased sample, which is critical for understanding the galaxy evolution and energy balance in 
a cosmological context. In this subsection, we discuss how the SF/AGN fraction varies with 
wavelengths, luminosity and redshift. 



4-3.1. SF/AGN fraction at mid- and total IR 



In Section 4, we have derived the 15 /xm luminosity function for the entire sample at 
z<0.3, as well as the 15 /im star-formation luminosity function. Integrating the differential 
LF, we estimate the luminosity density at 15 /im to be 1.5±0.3x10''LqMpc~^ to which 
the SF contribution is 9.0±2.2x 10^ Lq Mpc^^. This gives a SF fraction of ~58±19% to 
the integrated 15 /im luminosity density. Th en we extrapola te our results at 15 /im to the 
total IR. Following the technique described in lWu et al.l ( 120101 ). we convert 15 /im luminosity 
density to the total IR luminosity density for the AGN component and the SF component 
separately. Although the uncertainty in the SF fraction in total IR will be significantly higher 
because the error in converting Lis^^ to Lir also comes into play, this is still a critical quantity 
to obtain, especially for studies of the distant universe, where the PAH features are more 
difficult to measure, or the PAHs in high-redshift galaxies might have different properties (e.g. 
larger EWs) for similar dust mass fractions. We find that the star-formation contribution 
comes up to 83% of the total IR luminosity density. This is understandable as the SED of a 
star-forming galaxies is much steeper than AGN, thus the FIR emission will be dominated by 
SF. If we convert the IR luminosity for the star-formation component to the star-formation 
rate and integrate over cosmic time, we find our deri ved star formation r ate d ensity are 
consistent wi t h the dust extinct i on co rrected values of iHopkins fc BeacomI (120061 ) (see also 
Madau et al.l (119961 ): iLilly et al.l (Il996[ )). We repeat the same exercise at 24 /im. Integrating 
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the 24/im LF, we find the luminosity density at 24 /xm to be 3.4±l.lxlO^L0 Mpc~^ while 
the SF luminosity density is 2.6±l.lxl0^ Lq Mpc~^. This gives a SF fraction of ~78±42% at 
24 /im, higher than the number at 15 fim, while consistent with the concept that SF becomes 
more dominant at longer wavelength. We again convert the integrated 24 /xm luminosity 
to total IR luminosity density for star-forming systems and AGN separately and find the 
star-formation fraction to increase to ~89%, in agreement with the estimate from 15 /zm. 
If we take the average of the two (86%), then the SF fraction in the t otal IR luminosity 
density we have derived is higher than the study of IVeilleux et al.l (120091 ) for local ULIRGs, 



while more consistent with iNardini et al.l (120081 ). However, we need to bear in mind that 



both Veilleux and Nardini study the ULIRG population, whil e the lumino s ity of our sample 



is mostly in the range of 10 Loto 10 L0. Recent study of iPetric et al.l (120101 ) derived an 



average AGN fraction of 12% in the total IR for LIRGs, consistent with our estimate. 

As already noted (Figure [8]), the 15 /im SF LF departs from the 15 /xm LF most signifi- 
cantly at the bright end, an d this is again observed at 24 /xm. It su ggests that the SF fraction 
is a function of luminosity (lYuan et al.ll2010l : iHopkins et al.ll2010l ). In order to quantify how 
this fraction varies, we plot on the left panel of Figure M the SF fraction at 15 /xm versus 
the 15 /xm luminosity. We divide our sources into several luminosity bins and estimate the 
contribution of SF in each luminosity bin. The error bar represents the Poisson noise in 
each bin. We find that the SF fraction decreases as the 15 /im luminosity increases. This 
is expected since AGN plays a more important role in the energy budget for more distant, 
thus more luminous sources in our sample (see Figure [TU]) . On the right panel of Figure [HI 
we plot the SF fraction at 24 fim versus the 24 /im luminosity. We again observe a trend of 
the SF fraction decreasing with larger L24^m, however, the decline is much milder and the 
SF fraction is also higher at 24 /xm as compared to the 15 /xm. Finally, we convert the 24 /xm 
SF luminosity and 24 /xm luminosity to LgpiR and Lir, and show on the left panel of Figure 
[TUl how the total IR SF fraction varies as a function of IR luminosity, while on the right 
panel, we display the SF fraction versus redshift. We find that the SF fraction decreases with 
redshift, while there is little dependence of SF fraction with Ljr. The decrease of SF fraction 
with redshift can be understood since our sample selects a higher fraction of AGN-dominated 
sources as redshift increases (see also Figured]). We observe a mildly decreasing, or rather 
flat correlation between the IR SF fraction and Lir because: 1) At z<0.3, the luminosity of 
our sources only ranges from lO^Lo to lO^^L© and the majority of the sources included in 
this study are SF dominated; and 2) even for a source dominated by a powerful AGN in the 
mid-IR, its FIR emission could still be powered by star-formation, thus SF dominates the 
total IR luminosity. Because of the dominant contribution of SF in the FIR luminosity, we 
do not observe a strong dependence on Lir for the SF fraction. This suggest that the mid-IR 
might be a more reliable wavelength if one wants to study the SF/AGN fraction. The launch 
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of the Herschel Space Telescope has opened a new window for observi ng the cold dust i n 
the universe. Data from large area Herschel surveys, such as HerMES (lOliver et al.ll2010l ). 
will provide essential constraints on the FIR emission for sources in our sample, allowing a 
direct decomposition of SF/AGN in the total IR luminosity. Future studies using 70-500 /im 
data from Herschel will help to further constrain the SF fraction in the total IR luminosity 
density, and the associated uncertainties. 



4-3.2. Comparison with LF at z=0. 7 



In this section, we compare the LFs derived from the 5MUSES sample with relevant 
work at higher redshift to understand the evo l ution effects. The median redshift for our LF is 
0.12. When compared with 15 /im LF by IXul (120001 ) from ISO work, we find in general good 
agreement between our LF and Xu's LF, while at the high luminosity end, the two LFs show 
some discrepancy. Because of the higher median redshift of the objects in the high luminosity 
bin of our LF, this small discrepancy with the local LF hinted at evolution effects as a function 
of redshift. We then compare our work with studies at high redshift. 15 or 24/im LFs and 



LFs for star-forming galax i es have been deriyed by several differen t groups (ILe Floc'h et al. 



20051 : iMagneUi et al.ll2009l : iRuiopakarn et al.ll2010l : IFu et al.l 120101 ). In addition to deriving 
the rest-frame 15 /im LF from the 5MUSES sample, the availability of 5—35 /im IRS spectra 
also allowed us to decompose the star formation and AGN contributions to the IR spectral 
energy distribution in a more precise way, and thus derive the SF LF using the SF luminosity 
in each object from our sample. The most relevant work at higher redshift is at z=0.7 by 
Fu et al.l ( I2OIOI ). These authors use the IRS spectra of a z=0.7 sample to estimate the SF and 
AGN contribution to LF. In Figure [TTl we compare the 15 / im SF LF (red diamonds) from 
5MUSES with the corresponding 15 /^m SF LF at z=0.7 by IFu et al.l (120 lOh (blue squares). 
On the same figure, we also overplot the total 15 /im LF deriv ed from 5MUSES samp le (black 
filled circles) and its counterpart at z=0.7 from the work of iLe Floc'h et al.l ( 120051 ) (yellow 
crosses). We observe strong evolution effects in both the total and the SF LFs. With data 
only in two redshift bins (z=0.12 and z=0.7), we were not able to place stringent constraints 
on the amount of evolution on den s ity (a p) and luminosity (ol), however, if we adopt the 
values proposed by lLe Floc'h et al.l (120051 ) and evolve our total 15 /im LF at z=0.12 to z=0.7 
by a factor of aD=2.1 and aL=2.6, we find a good match between the evolved LF (black 
dotted line) and the observed data (yellow crosses). 

We follow a similar approach to evolve the 15/im SELF from z=0.12 to z=0.7. In 
Figure [T2| we show that if we evolve our 15/im SFL F at z=0.12 (black solid line) by a 
factor of aD=2.5 and q;l=2.6, it matches very well (reduced = 0.6) with observed rest- 
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frame 15 /iin SF LF at z=0.7 from lFu et al.l ( l2010l )'s work (blue dashed line). Could a pure 



density evolution or luminosity evolution explain the difference we observe in the SFLF at 
z=0.12 and z=0.7? In Figure [T^ we show the best fit to the SF LF at z=0.7 if we only allow 
density evolution on our local SFLF (green dash-dot-dot line). The best fit returns aD=5.3 
and the reduced x^=65.6. On the same figure, we also show the best fit when only luminosity 
evolution is allowed (yellow dotted line). We find aL=4.3 and the reduced x^=70.0. Even 
though we only have data in two redshift bins, which makes it difficult to place stringent 
constraints on the amount of evolution we require on luminosity and density, we can at least 
confirm from our work that neither pure luminosity nor pure density evolution is sufficient 
to explain the difference between SF LFs at z=0.12 and z=0.7. 



We should however point out some caveats. Although both our work and iFu et al. 



( l2010l )'s work used IRS spectra as diagnostic tools for distinguishing SF and AGN, our LF 
is derived by separating the energy density contribution in each source, while they use IRS 
spectra to determine the major source of energy in each object, and then remove the AGN 
dominated sources when building their SFLF. Because they only have IRS spectra for ~ 40 
galaxies from their flux limited sample, the low end of their SFLF is obtained by shifting the 
MIPS 24 /im luminosity to rest-frame 15 /im (z=0.7), and they assign these low luminosity 
objects as dominated by star formation. Given our very rigorous analysis, we clearly see 
some differences in the total LF and the SFLF at 15 fim even at the low lur ninosity e n d (see 



the comparison of the black dot-dashed line and the black solid line), thus iFu et al.l (j2010[ ) 
might have overestimated the SFLF below Li5^m=lO^°'^L0. 



5. Conclusions 

We have analyzed a Spitzer spectroscopic sample of 24 /xm selected objects in the SWIRE 
and XFLS fields to derive the 15 and 24 /im luminosity functions at 0<z<0.3. When com- 
bined with local and high-redshift studies, this provides critical information on understanding 
the evolution of energy budgets over the past three billion years (z<0.3). Our conclusions 
are summarized as follows: 

1. We have derived 24 and 15 /im LFs for our sample at z<0.3. The availability of the 
5-35 /im allows us to make K-corrections directly using the observed SED. The 24 and 15 /tm 
LFs display rather shallow slopes at the bright end, which is due to the increase of AGN 
contribution in more luminous systems. 

2. Using the 5-35 /im IRS spectra, we have decomposed the 5MUSES objects into star- 
formation and AGN components. The SF fraction is taken to be the median likelihood of 
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the probability density function. We calculate the 15 and 24 //m SF luminosities in each 
object, and subsequently build the SF LFs. The SF LFs can be described with a Schechter 
function. 

3. We have also estimated the SF contribution to the integrated 15 and 24 //m luminosity 
density for our sample. The SF fraction is found to be 58% at 15 //m and it goes up to 78% 
at 24 nm. Using the conversion factor from Lmir to Lir for star-forming galaxies and AGNs 
respectively, we found the SF fraction to be ~86% in the total IR luminosity density. 

4. The SF fraction is also found to be a function of luminosity/redshift, decreasing as 
luminosity or redshift increase, while the trend is more obvious in mid-IR, suggesting that 
mid-IR wavelength is more sensitive to the presence of AGN. 

5. Both luminosity and density evolution are required to explain the difference in the 
observed SFLF between this sample and similar studies at z~0.7. 
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Table 1. mid-IR Luminosity Function Derived from the 1/Knax Method 
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Fig. 1. — The Redshift distribution for 309 out of 330 sources in the 5MUSES sample (dotted 
hne). The blue, yellow and red lines represent the distribution for SB, composite and AGN- 
dominated sources. The dashed vertical line indicates the redshift cut of z<0.3 on which 
this paper is focused. 



Table 2. Results of the FITs to the mid-IR Luminosity functions 
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Fig. 2. — The ratio of flux densities at SDSS r band to the MIPS 24 /xm band versus 
redshift for 5MUSES sources. The blue circles, yellow diamonds and red crosses repre- 
sent the SB, composite and AGN sources. The green stars represent sources without red- 
shifts. The dotted vertical line indicates our redshift cut of z<0.3. The dashed line indicates 
log[f^(r)/f^(24^m)]=-2.6. All 13 objects with log[f^(r)/f^(24/im)]<-2.6 are located at z>0.3. 
17 out of 21 unknown-z sources have log[fi,(r)/fi,(24/im)]<-2.6, and are thus most likely to 
lie at z>0.3. 
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Fig. 3. — The completeness correction factor for 5MUSES versus redshift. Instead of using 
an average correction factor of ~3.4 (1111/330), we benchmark our sample with the XFLS, 
where redshifts are available for this field. We first derive the fractional contribution of 
number counts in different redshift bins for the XLFS field at Ij, >5mJy. Using this as 
a reference, we predict the number counts in each redshift bin for the 5MUSES sample. 
Then we divide the predicted number of objects by the observed number and obtain the 
completeness correction in each redshift bin. Finally, we use a second order polynomial to fit 
the data and obtain the completeness correction factor as a function of redshift. The error 
bars represent the Poisson noise. 
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Fig. 4. — a) Upper panels: Examples of our decomposition on the IRS spectra for typical 
starburst, composite and AGN-dominated sources in 5MUSES. The black solid line is the 
IRS spectrum and the green shaded region indicates the uncertainties. The blue dotted line 
is the scaled star-formation galaxy template; the yellow dashed line is the scaled quasar 
template; the red solid line is the best-fit composite spectrum. The 6.2 fim PAH EWs (a 
negative value indicates an upper limit), redshifts of the objects, as well as the reduced 
values are also shown on the plot, b) Lower panels: The probability density function for the 
star-formation fraction at 15 /xm for each corresponding object. The median fraction and its 
16th to 84th percentile range is also indicated. 
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Fig. 5. — Left panel: The SF fraction at 15 jj,m versus the 6.2 jj,m PAH EW. The dotted 
Hne indicates 6.2 /xm PAH EW=0.2 fim. Middle panel: The SF fraction versus the mid- 
IR continuum slope of fi,(24/im)/fi,(15/im). Both the PAH EW and continuum slope are 
correlated with the SF fraction at 15 //m with scatter. The two sources on the middle panel 
with very low f,^(24//m)/f,^(15//m) ratios and high SF fractions have their mid-IR spectra 
dominated by PAH emission, while their continuum slopes do not rise quickly. This also 
indicates that continuum slope alone has a high uncertainty when it is used as a SF indicator. 
Right panel: The SF fraction versus the infrared luminosity. There is no clear correlation 
between these two parameters. Note that at the high luminosity end, we don't observe many 
objects with low SF fraction. This is because we have not included very high luminosity 
sources in this study (z<0.3). 
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Fig. 6. — The rest-frame MIPS 24 /im luminosity function at z<0.3 derived from the 5MUSES 
sample using the 1/Vmax method (black circles). The yellow squares represent the LF if we 
were to apply a uniform completeness correction factor of ~3.4. The blue di amonds and green 



crosse s represent the 24 /xm LFs at 0.05<z<0.2 and 0.2<z<0.35 derived by lRujopakarn et al. 



(I2OIOI ). The dashed line is the fit to the data points of our 24 /xm LF assuming a double 
power-law exponential function. The knee (L*) of the LF is represented by the red star on 
the plot. The dotted line denotes the luminosity calculated at the median redshift of this 
sample corresponding to the flux limit at 24 yum. 
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Fig. 7. — Top panel: The monochromatic 15 /im luminosity function at z<0.3 derived from 
the 5MUSES sample (filled circles). The dashed line is the fit to the data points assuming 
double power-law exponential function. The blue dot-dashed line represents the fit with a 
Schechter function to our data. The knee (L*) of the LF is represented by the red star. The 
dotted line corresponds to the luminosity calculated at the median redshift of our sample for 
a galaxy at the flux limit. The local 15 /im luminosity function from IXul (120001 ) is overplotted 
with the yellow diamond for cor npar i son. Our 15 /xm LF is in excellent agreement at the faint 
end with the ISO 15 /xm LF by IXul (|2000[ ). but is slightly higher at the bright end. Bottom 
panel: The median redshift in each luminosity bin we have used to derive our LF. The solid 
line represents the median redshift of Xu's sample. 
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Fig. 8. — a) Left panel: A comparison of the 15 fim luminosity function (black filled circles) 
and the 15 /xm star-formation luminosity function (red squares). The dashed line and the 
dotted line are the f its to the two LFs. The blue crosses represent the local AGN-excluded 
LP by iPozzi et al.l ( 12004 ) and appear to be in good agreement with our 15 fim SF LP. 
The green diamonds represent the AGN LF derived from our sample and the dot-dashed 
line provides a fit to the AGN LF with a double power-law exponential function. b) Right 
panel: A comparison of the MIPS 24 /xm luminosity function (black filled circles) and star- 
formation luminosity function (red squares). The dashed lines are the fits to the total LFs 
with a double power law exponential function, while the dotted lines are the fits to the SF 
LFs with a Schechter function. The blue diamonds and green cros ses represent the 24 /im 
SF LFs at 0.05<z<0.2 and 0.2<z<0.35 bv iRuiopakarn et al.l ^01^ . 
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Fig. 9. — a) Left panel: The fractional contribution of the star-formation luminosity in each 
luminosity bin at 15 /im. The dotted line is a second order polynomial fit to the data, b) 
Same as a), but for 24 //m. 
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Fig. 10. — a) Left panel: The fractional contribution of the star-formation luminosity in 
each luminosity bin in the total IR luminosity, converted from the results at 24 /xm. The 
dotted line is a second order polynomial fit to the data, b) The SF fraction in the total IR 
luminosity versus the redshift. The dotted line is a second order polynomial fit to the data. 
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Fig. 11. — Comparison of LF at < z >= 0.12 from our sample with LFs at < z >= 0.7. 
The black filled circles represent the total 15 fim LF and the red open diamonds represent 
the 15 fim SF LF derived from our sample. The black dash-dotted line and the solid line 
are the best-fits to the SFLF and total LF derived from the 5MUSES s ampl e at 15 /xm. The 
blue open squares represent the 15 fim SF LF at z=0.7 fromlFu et al.l (120101 ) and the yellow 
crosses represent the total 15 /xm LF from iLe Floc'h et al.l (120051 ). The blue dashed line is 
the best-fit Schechter function from iFu et al.l (120101 ) to the 15 /im SF LF at z=0.7. The 
black dotted line is derived by e volving the luminosity and density of LFs from 5MUSES by 
a factor of ai^=2.6 and ar.=2.1 (ILe Floc'h et al.ll2005r) and it appears to match well with the 
15 ^m LF at z=0.7 bv lLe Floc'h et al.l J2OO5I 1 
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Fig. 12. — - Compar i son o f the 15 /iin SF LF from 5MUSES (red diamonds) aX < z >= 0.12 
and from iFu et al.l (120101 ) (blue squares) aX < z >= 0.7. The black solid line and the blue 
dash-dotted lines are the best fits to the two LFs respectively. The green dash-dot-dot line 
represent the best fit when only evolution on density is allowed (aD=5.3). The yellow dotted 
line represents the best fit when only evolution on luminosity is allowed (aL=4.3). The black 
dashed line represents the best fit when both density and luminosity evolution is allowed 
(aD=2.5, aL=2.6). 



